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Introduction 


The Ukrainian Carpathians is the territory with a great threat of floods. This is due to natural and climatic 
conditions of this region characterized by mountainous terrain, high density of hydrological network and a 
significant amount of precipitation. Amount of precipitation here ranges from 600 mm on plains to 1600 
mm on mountain tops. The main factors of flood occurrence are excessive amount of precipitation, low 
water permeability of soil and a high proportion of low-permeability rocks (flysch layers with a 
predominance of clay layers). Therefore, catastrophic floods in the region were also observed in previous 
centuries, when the anthropogenic impact on the environment, including forest ecosystems, was not 
comparable with the current one (Figure 1) (Mandryk et al., 2017; Arkhypova et al., 2019). 





Figure 1: Flood in Ivano-Frankivsk region (June 2020). Courtesy: Photo from open sources 


The rivers of the Carpathian which have flood hazard area belong to transboundary waters. In the basins of 
Carpathian rivers, some part of which belong to neighboring countries, the European Union Water and Flood 
Directives are implemented in accordance with the Association Agreement between Ukraine and the 
European Union. These directives are implemented with the funding from European Union (Prykhodko et 
al., 2020; Kinash et al., 2019). 


The estimation of the level of flood waters is an urgent scientific and technical problem in the context of 
complex geo-climatic processes, which lead to catastrophic floods, taking place in European countries (such 
as the floods of Ukraine in 1927, 1941, 1969, 1974, 1980, 1998, 2001, 2008, 2010 and 2020, of Poland and 
Slovakia in 2012, of the Czech Republic and Great Britain in 2013, and of Germany in 2002). The stated 
problem is addressed by a wide range of scientists studying and solving the applied problems by using 
precise models of filtrating flows (Oliynyk and Panchuk, 1992; Oliynyk and Steyer, 2012) and empirical 
dependences (Mandryk, Pukish and Zelmanovych, 2017; Maslova and Susidko, 2006; Sosedko, 1980; 
Leontiev, 2009), which are established on the basis of experimental data analysis (Zasidko et al.,2019; 
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Pietrzak et al., 2018). However, the above-mentioned works do not suggest the idea of the possibility to 
build the express method for forecasting the phenomenon of flood development. 


In the present paper, based on the use of methods for approximating the curves and the results of 
experimental studies of flood waters, a method of mathematical description and forecasting of the flood 
development is suggested. Presented method gives next possibility: 


e to apply the mathematical methods for studying of real processes (floods, seasonal phenomena of 
emergency water distribution; 


e to offer the functional analytical structure based on the statistical information about the floods in 
concrete region, which give the possibility to predict the flood’s duration and intensity; 


e torestore the ordinary differential equation which describe the flood process development using the 
information about one’s solution; 


e to choose the factors that have the significant influence on the flood’s characteristics; and 
Methodology 
Mathematical modeling of the flood development process 


When analyzing the real phenomena of flood situations, it is possible to depict schematically the relationship 
between time (1 ) and the level of flood waters (/ ) (Figure 2). 











Figure 2: Diagram of the relationship between time and the level of flood waters 


Any flood is characterized by a period of development a), a period of its critical (maximum) intensity b) 
and a period of decline c). The scheme shown in figure 2 does not consider cases when for a short period of 
time the flood has several peaks. In this case the phenomenon shown in Figure 2 can be simulated by 
application of several floods with one peak. It 1s necessary to reproduce the analytical structure of the 
function shown in Figure 2. From a mathematical point of view, this function must be characterized by the 
following conditions: 


f (0) =0; 
Ixy: f (0)=0; f (%) <0; (1) 
lim f(x) =0. 


Obviously, many functions meet such conditions, and, therefore, the problem of constructing the flood 
development model from a mathematical point of view 1s ill-defined (Leontiev, 2009) to regulate. It is 


Oleg Mandryk, Andriy Oliynyk, Roman Mykhailyuk, Lidiia Feshanych 67 


Grassroots Journal of Natural Resources, Vol.4, No.2 (June 2021), pp.65-76 |ISSN 2581-6853 | CODEN GJNRA9 
Doi: https://do1.org/10.33002/nr258 1.6853.040205 OPEN GACCESS Freely available online 


needed to develop certain algorithms using additional information about type of the function that would 
satisfy the conditions (1) (Tikhonov and Arsenin, 1979; Zorich, 1981; Filippov, 2007). 


The simplest of the known functions is the following one: 
(Mate, (2) 
the graph of which is shown in Figure 3. 
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Figure 3: The graph of the function f(t) =t-e" 


Verified by means of direct check, it has been established that, for the function (2), all three following 
conditions (1) are met: 


L. f(0) =0-e° =0: 

2 AF) if) = 0: 
(t-e )=e —-e'+t-e' =O>5treHl; 
f @Q=—-e'-e't+t-e* =(-2412)-e"; 
f d=-e' <0; 


. . 1 
3. lim — = lim— =0. 
to go! t>0 g! 


However, application of function (2) to describe flood phenomena is connected with some problems. It 
cannot be changed, and its type cannot be selected according to certain results of experimental studies and 
analysis of statistical data on floods. To solve this problem, the following function can be suggested: 


y=f(=t-e%, a>O, (3) 
which obviously also satisfies the stated conditions (Figure 4): 
1. AsO "20% 
2. a7) a0: 
(t-e“) =e “—a-t-e“ =(l-a-t)-e” =0; 


t=—; 
a 
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f Q®=-a-e“+(l—a-t)-(-a)-e“% =e“: (-2-a+a’-t)=0; 


" i 
f (=-e':(-2-a+a)=-ae' <0; 
a 


=0. 





a. lina =i : 


t— oo eo t— oo a: eo 
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Figure 4: The graph of the function y= f(t) =t-e“, a>O 


Dependence (3) has a parameter a >0 that allows to obtain a whole range of curves, which, however, are 
topologically (according to their spatial location) similar, and that does not allow to introduce many 
parameters on which the level of flood water rising depends. That is why, the following two-parameter 
model is proposed to determine the function of type (1) of the form (Figure 5). 


y=t"-e“,a>0,n>0, (4) 
for which all conditions (1) are determined: 
et a : : a oe n 
1. (t"-e%) =n-t"'-e%-a-t"-e% =t" | -(n-e“+a-t-e“)=t" -e"-(n-a-t)h=05t=—-; 
a 
(t"-e“) =-a-e“-t"' -(n—-a-t)t+e-t” * -(n—-1)-(n—a-t)+e-t" | -(-a) = 
=e 7"? [(-a)-t-(n-a-t)+(n-1)-(n-a-t)-a-t]=e- 1"? | -a-t-nta’-t? |+ 
+(n —n-n-a-tta-t)—a-t]=e"-1"?-|-a-t-nta’-t? +n -n-n-a-t |= 


=e "fp" | a’ tt —Q-atntn -nj=e*r"? | (a-t-ny -n | 


f Q=-e")"? [=n -n]=-n-e "Oy"? <0; 
a a a 





2. f(0)=0; 
a. lim a =0 according to the L'Hospital's Rule (Zasidko et al., 2019). 
[>0 @ 
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Figure 5: Family of graphs for the function y=?"-e“, a>0,n>0 


By applying, for example, dependence (4), it is possible to restore the differential equation of the process 
described by each of the functions (2) — (4). 


To fulfill this, we can use the following property of linear systems: if the function y,(x) satisfies the first- 
order differential equation of the form f(y,y)=0, then the following condition is met (Berestneva, 
Marukhina and Shevelyoev, 2012): 

yx) y (x) 


_ [=0. (5) 
yx) y,() 








Then, by substituting in the determinant (5) of the function y,(x), according to the properties of 
determinants, the equation (5) is satisfied identically. Applying similar values, it is possible to obtain 
differential equations, the solutions of which are the corresponding functions: for the function (2), equation 
(5) requires the form: 








y: y _() 
e'-t-e' t-e' 
or after summing up the similar addends: 
y t-(l-t)-y=0, (6) 
for the function of type (3), the form: 
y Dy. _0 
t : ee ae a t “a: e* 
or: 
y -t-(l-t)-y=0, (7) 
and for the function of type (4), the form: 
y y = 
t” po n- t” | —at af" -d as 








where the following values are obtained: 
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Y gilt —a n —a nN 
y-t"-e" —-t"-e C-a)-y=0; 


y =(2-a)-y. (8) 
ft 


For each of the equations (6) — (8), especially for equation (8), there can be set the initial conditions of the 
form: 


W(t) = Yo» (9) 
setting the conditions in the form of (9) allows to take into account the initial level of flooding: 


—a az —n n ft = _ 
y=y-e 0 Le a? -t -»(<] oo to) 


It is also important to establish the physical content of coefficients n and a, which may be functions of the 


form: 
‘ =n(x, N- bes Ms t); 


a =a(x, are t), 


(10) 


where the values are soil permeability, air humidity, water intake humidity, relief features, etc., and the 
variable is time. 


The question arises how to determine the numerical values of parameters (10) for different types of floods 
and how to choose the variables related to (10), which most significantly influence these dependencies. It 1s 
also necessary to define the components of the equation (8). Function y(t) represents the level of flood 
waters at some point of time, which varies in proportion to itself with an alternating-sign coefficient of the 
form: 


k=" -a, (11) 


where n and a are empirically determined coefficients. When k >0, the intensity of the flood increases; 
when k <0, it decreases. 


Let the statistical data on its course (y,, t,) are known for the chosen flood. According to these values 
(y,,t,) i=L..,N, where Nis the number of observations, which during intense floods can be quite 


significant due to the fact that flood water levels are monitored regularly, estimation of parameters n and 
a can be done by using a linear regression apparatus. For this purpose, function (4) can be written (applying 


logarithm problem) as: 
In y=nInt-—at, 


Iny Int (12) 
— =n—- 4, 
t t 


Int 
introducing the notation = eas fle- es kn; Y= —a, we obtain: 
t t 
Go= Khe bf, (13) 


that is, the linear regression equation. Using the known formulas (Zamikhovskii et al., 2014) for linear 
regression coefficients, we obtain y -t—(1-t)-y=0: 


fo N >. Hle-) | Hoy fe 
NY - (Sy (14) 


, | 
dee Oe — RS ff, 
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from where, using communication formulas, is obtained the following: 

n= ke f= —a, (15) 
that is, parameters n and a can be determined unambiguously, that allows us to talk about the constructed 
regulable algorithm for the ill-defined problem of restoring the function under conditions (1). To determine 
the factors influencing the process, the method of (Pietrzak et al., 2018) associative analysis is used, which 
allows to identify variables and parameters that affect the process of increasing flood waters. Suppose, 
according to the results of experimental research and analysis of statistical data, a quantitative characteristic 
of some parameter x, has been established, the range of change of which can be divided into two segments 


that correspond to approximately equally probable values x,. Simultaneously, the following table is 
arranged: 





where f stands for the value of the level of flood waters, f, stands for some average value that divides the 
range of flood waters change into intervals in which the values f are distributed approximately equally in 
number. A, B, C, D are the numbers of comparison results that correspond to the specified values f and x, . 


The following values are calculated: 
A+B=n;; 
C+D=n,; 
A+C=n,; 
B+D=n,. 


Obviously, the total number of experiments is equal to either n,+n, or n,+n,. The contingency ratio is 


calculated according to the formulas: 


AD— BC 


If w > 0,3, then the relationship between the values is considered to be proved and significant; it should be 
studied in more detail, but if yw < 0,3, then the relationship between the specified values can be considered 
insignificant. 


Results 
Methods of forecasting 


The following method of estimating the level of flood waters and forecasting their development has been 
proposed. Suppose that by studying the floods that took place in the region under study with the help of 
methods (12) — (15) the ratio: 

Pa ee ea aks (17) 
where k is the number of the studied floods that has been obtained. 


The parameters x,,x,,...,x, , that affect the formation and development of floods, are determined. They are 
water and physical properties of soils (water permeability, water intensity), air humidity, wind direction and 
speed, etc. It is considered that based on experimental studies, these values are known for each of the floods 

x',,X',,..x'. According to 


m e. 
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the method (Pietrzak et al., 2018) of associative analysis (16), for each of the values x,, s=1,...,.M the 
contingency ratio (16) 1s calculated and the level of connection between the corresponding x, and f is 
determined. Thus, the number of values x,, i =1,...,M that influence f , 1s reduced, and in the future only those 
x,, j=1,...,M,), M, <M, that influence the process are considered. The presented results allow to optimize 
the experimental research procedure — the number of parameters x et = eee ys for which it is necessary to 


develop methods of experimental evaluation and control, is reduced. When studying the possibility of 
flooding in the given region, the values x’,, j =1,...,.M, are determined; after that the formula (17) and the 


parameters x,, j=1,...,M, determined for the floods f, that influence the floods f,, are found out. 


The following value 

argmin{ >) (x) xP) f= JF, (18) 
is determined in order to obtain more accurate forecast, the values j,,s =1,....N,,N)(N, = 2) that are the 
nearest to j are determined. We choose those formulas (17) that correspond to the determined j, . Applying 
the corresponding dependences f, (17) determined for j, (in the simplest case j, is the only one and 


j, =j , determined based on (18)), along with corresponding graphs, we choose the possible level of 
flooding and its intensity. 


The scheme for estimating the level and duration of floods is estimated, for example, as shown in the graph 
(Figure 6). The values and “are determined on the basis of statistical analysis. 


iNT er ik ets ms et 





Figure 6: Scheme for estimating the flood level and duration [Explanation to figure 6: y, stands for critical 
increase in the level of flood waters, At for duration of a flood. ] 


Conclusion 


The method of estimating the flood water level and forecasting of flood phenomena has been suggested. 
Compared with the existing methods, the presented one is the method allows to offer the function analytical 
structure based on the statistical information about the floods in concrete region, which give the possibility 
to predict the flood’s duration and intensity, to restore the ordinary differential equation which describe the 
flood process development using the information about one’s solution. It 1s possible to change the order of 
differential equation to describe the processes more in detail, to apply the method of associative analysis for 
the choosing the factors which have the significant influence on the flood’s characteristics. From the one 
side it is the method with the strict mathematical justification but from the second side it is quite simply in 
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realization in correspondent institution which activity is devoted to the struggle with the floods in our region. 
The direction of further research may be related to the development of experimental means to determine the 
parameters that affect the process of flood formation. 
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